Gapless interface states at the junction between two topological insulators 
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We consider a junction between two topological insulators, and calculate the properties of the 
interface states with an effective low energy Hamiltonian for topological insulators with a single 
cone on the surface. This system bears a close resemblance to bilayer graphene, as both result from 
the hybridization of Dirac cones. We find gapless interface states not only when the helicity direction 
of the topological surface states are oppositely oriented, but they can also exist if they are equally 
oriented. Furthermore, we find that the existence of the interface states can be understood from 
the closing of the bulk gap when the helicity changes orientation. Recently, superluminal tachyonic 
excitations were also claimed to exist at the interface between topological insulators. However, here 
we show that these interface states do not exist. 
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I. INTRODUCTION 

Topological insulators (TIs) 1-12 are a newly discovered 
class of materials, that have attracted a lot of interest 
from the condensed matter community in the last few 



years 



13-18 



The strong three-dimensional TI, e.g. Bi 2 Se 3 , 
is insulating in the bulk with gapless surface states pro- 
tected by time-reversal (TR) symmetry as a consequence 
of band inversion by strong spin-orbit coupling (SOC). 
Elastic backscattcring of these states is forbidden due to 
TR symmetry, and they remain gapless for any TR in- 
variant perturbation of the system that does not close 
the gap. At low energies, the surface states of Bi2Se3 
are given by a single Dirac cone, i.e. they have linear 
dispersion and a helical spin texture 9,11 . 

Recently it was shown that there should also exist pro- 
tected gapless interface states, based on symmetry argu- 
ments, at the junction between two TIs with opposite 
helicity direction 19 . These interface states are not pro- 
tected in the same manner as the surface states of a single 
TI, instead they are protected by mirror symmetry. The 
physics of this system resembles certain aspects of bilayer 
graphene, because both result from the hybridization of 
Dirac cones. 

In this paper, we make a systematic study of the pos- 
sible combinations of TIs using a quantitative model for 
a strong TI 9,20 , and we show the existence of different 
types of interface states. This model can be derived from 
fc • p perturbation theory in which the full Hamiltonian is 
projected on the subspace of states that dominate near 
the r point. We follow the same approach as a very 
recent paper 21 that claims the existence of tachyonlikc 
interface states in this system. We found, however, that 
this model does not predict tachyonlike solutions. The 
reason is technical, and due to a wrong implementation 
of the model under certain pathological conditions. 

The paper is further organized as follows. In Sec. II we 
briefly discuss the model that is used to describe the TIs, 
and we explain how the interface states were calculated. 
Next, we present and discuss our results in Sec. Ill, and 
we prove that the tachyonlike interface states are not 
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FIG. 1: (color online). Junction between Til and TI2 with 
opposite helicity direction. The spin of the surface Dirac cone 
lies in the plane, perpendicular to the direction of propagation 
(given here by the y-axis), oriented along the positive (green) 
and negative (red) 2-axis for Til and TI2 respectively. 



physical solutions. 



II. MODEL 



General 



The low energy physics of a bulk strong TI, like 
Bi2Sc3, is described up to order fc 2 by the effective 
Hamiltonian 9,20 
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where cr are the Pauli matrices for the Bi and Se sites, 
and 

£(fc,fc z ) = C + Cifc 2 + C 2 A: 2 , 
M(k, k z ) = M Q + M x k\ + M 2 k 2 , 
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(k x ,k y ) and k± 



k x i ikn 



The 



Schrddinger equation is given by HQ = where 
<I> = (Bi-f, Sef, B14., Sej,) gives the amplitude of the spin- 
orbit coupled p z orbitals that mainly contribute to the 
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low energy physics of Bi2Se3. The band inversion which 
characterizes the topological phase, is determined by the 
condition M Mi < (Ref. 20). 

Surface states obtained from this model have linear 
dispersion near the T point with Fermi velocity v = 

1/2 

{\A\/h) [l - (Ci/Mi) 2 ] 1 , and a helical spin texture 22 . 
The hclicity direction of the spin texture depends on the 
relative sign of the parameters A and B, i.e. (s x . y ) ~ 
sgn (AB) k Vt x (Ref. 20). 



B. Junction 

We consider a junction at z = between two TIs, 
Til (z < 0) and TI2 (z > 0) (Fig. 1). The junc- 
tion breaks translation symmetry in the z-direction, and 
we let k z — !• — id z to obtain a system of second order 
homogeneous differential equations H(k, —id z )Q = EQ 
with $ = e lk ' p 4>(z). This is solved with the ansatz 
4>(z) ~ e Xz, $(E, k), where ^(E,k) is an eigenvector of 
the Hamiltonian. We obtain a system of algebraic equa- 
tions that has a nonzero solution if \H(k, —iX) — E\ = 0. 
Details on the solution method are given in the Ap- 
pendix. Solving for A yields four doubly degenerate 
X a {E,k) (a = 1,...,4) in general, each corresponding 
to two eigenvectors * SjQ (s = 1,2) (Eq. (A.2)). The to- 
tal solution is then given by the linear combination 



(2) 



where the coefficients C StCC (E,k) are found from the 
boundary conditions. 

We want to study states localized at the junction, and 
therefore we only consider the two A' 1 "- 1 with Re A' 1 - 1 > 
for Til and Re A^ 2 ' < for TI2, to construct the wave 
function (2) (sec Eq. (A.l)). Here m labels Til and 
TI2 respectively. The remaining boundary conditions are 
given by the continuity of the wave function and the cur- 
rent at the junction 23 
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k z — — idz. 



k z ——id z 



(3) 



The nonzero solutions C^a (E,k) ^ of this system of 
homogeneous equations, define the dispersion and wave 
functions of the interface states. We numerically solved 
this system on a (E, k) grid by means of the condition 
number and singular value decomposition of the coeffi- 
cient matrix. 



III. RESULTS 

The parameters of the model (1) are taken from Ref. 9 
and were obtained from ab initio calculations of Bi2Se3. 
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FIG. 2: (color online), (a) Dispersion of the interface states 
for B< 2 > = B«, and A< 2 > = -A« (1), = -2A« (2) and 
A (2) = -10A (1) (3). (b) Decay length of the interface states 
for the two bands from (a) with A^ 2 ' = —A' 1 ', (c) Density of 
the states of band lb from (a) marked with a cross, near the 
crossing point (solid), and closer to the bulk region (dashed). 



We want to study interface states between TIs with op- 
posite helicity, and therefore we focus only on the param- 
eters A^ 2 ) and B^ 2 \ which determine the helicity direc- 
tion of the surface states of TI2. We find gapless interface 
states if the helicity direction of the surface states of the 
TIs is equally as well as oppositely oriented. Specifically, 
we find solutions if the sign of the parameters A or B is 
opposite for the TIs (case 1 and 2), corresponding to TIs 
with opposite helicity, and also if the sign of both param- 
eters is opposite (case 3), corresponding to TIs with equal 
hclicity. We find no solutions for the other case in which 
both TIs have equal helicity, i.e. if both parameters have 
the same sign. To simplify the calculation and because 
the Hamiltonian (1) has full rotation symmetry about 
the z-axis, we take the wavevector along the x-direction. 



3 



*2 



1 

X 



A 

... 









0.1 0.2 0.3 

Wavevector (nnf 1 ) 



(2) f 




3-1.5 


V 

1.5 C 



> 

c 
c 

a 



0.3 



0.2 



0.1 



-0.1 



-0.2 



•1 



A 

... 









2 




0.1 0.2 0.3 

Wavevector (nnf 1 ) 



-1.5 1.5 

z (nm) 



FIG. 3: Dispersion of the interface states between TIs for 
B (2) _ _ B m and A (2) _ WA m The density ( arb units ) is 

shown on the right for the two states that are marked with a 
cross on the dispersion. 



FIG. 4: Dispersion of the interface states between TIs for 
B (2) = _ B 0) and A (2) _ _ 10y 4(i) The density (arb. units) 
is shown on the right for the two states that are marked with 
a cross on the dispersion. 



A. Case 1 

First wc consider the case where = and 

y 4( 1 )^4(2) < o. This corresponds with TIs that have sur- 
face states with opposite helicity direction (Fig. 1), and 
the magnitude of A^ 1,2 ^ determines the Fermi velocity of 
the surface states. In Fig. 2(a) the dispersion of the in- 
terface states is shown for three different values of A^ 2 '. 
The interface states have linear dispersion (in the radial 
direction) around k ^ which results from the hybridiza- 
tion of the Dirac cones from the separated TIs, similar 
to A A stacked bilayer graphene 24 . The cones intersect 
closer to the center if their Fermi velocity is smaller. 
From the energy scale in Fig. 2(a) one can see that the 
conduction/ valence band of the cone is pushed into the 
bulk region of the band structure, which has a band gap 
around 0.3 eV 9,11 . Correspondingly, as the bands ap- 
proach k = 0, the decay length = max(l/|Rc A Q |) di- 
verges exponentially, as shown in Fig. 2(b). 

The total density near the junction is shown in Fig. 
2(c) for a state close to the crossing point and a state close 
to the bulk region. Both states are marked with a cross 
in Fig. 2(a). We see that the character of the individual 
surface states is completely lost and the interface states 
are spread over the entire junction. 



B. Case 2 

Next we consider the case where B^B^ < and 
J 4( 1 ) J 4(2) > o. Again this corresponds with TIs that have 
surface states with opposite helicity direction (Fig.l), and 
the magnitude of A^ 1,2 ^ determines the Fermi velocity of 
the surface states. In Fig. 3 the dispersion of the interface 
states is shown for B^ = -flW and A^ = 10A (1 K Un- 
like case 1, there is little interaction between the surfaces, 



and the cones lie on top of each other. 

The density is also shown for two states in Fig. 3, one 
on each cone. We see that the states on each cone are not 
localized at the junction, but instead they are localized 
inside one of the TIs, as we would expect if there is little 
interaction, and the individual surface states keep most 
of their original character. TI2 has the largest Fermi 
velocity, correspondingly the states on the steepest cone 
arc localized in TI2. The Fermi velocity of the cones 
in the interface spectrum is two times as large as the 
original cone of the separated TI for Til and about 0.8 
times smaller for TI2, and thus the cones are attracted 
to each other. 



C. Case 3 

Finally we consider the case where B^B^ < and 
A^'A^ < 0. Unlike the previous two cases, this cor- 
responds with TIs that have surface states with equal 
helicity direction. The dispersion of the interface states 
is shown in Fig. 4 for B^ = -B^ and A^ = -lOA^. 
We see that we get a combination of the previous two 
cases, a linear spectrum near the center and a crossing 
away from the center. The density is also shown for two 
states in Fig. 4, localized on Til, one near the center and 
one near the crossing point. We see that near the center, 
we have the same behavior as in case 2, and closer to 
the crossing point the interaction between the surfaces is 
stronger as indicated by the shift of the density towards 
the junction. The states on the steep cone do not dif- 
fer much from those of the steep cone of Fig. 3 and are 
localized within TI2. 



FIG. 5: Schematic phase diagram as function of A and B. 
When either of the parameters changes sign, the gap closes at 
(k,k z ) shown and the helicity changes orientation. Crossing 
from region to region 1-3 corresponds with case 1-3. 

D. Discussion 

If the helicity direction of the TIs is opposite, the ex- 
istence of the gapless interface states can be understood 
from scattering of surface states at the junction between 
the two TIs, as shown in Refs. 19 and 25. At normal 
incidence, the incoming state on the surface of Til can- 
not scatter backwards because of spin conservation. This 
is similar to graphene, where pseudospin conservation re- 
sults in Klein tunneling through potential barriers at nor- 
mal incidence 26 . However, when TI2 has opposite helic- 
ity, the incoming state cannot go forward because there 
are no states in TI2 that conserve both spin and momen- 
tum. This paradox is solved by the existence of gapless 
interface states between the TIs. At normal incidence the 
incoming state has to go into the interface. Consequently, 
these states should exist for k x = or k y = 0. 

The origin of the surface states in TIs can be under- 
stood from the closure of the bulk band gap if the Hamil- 
tonian is transformed continuously from a trivial insula- 
tor to a TI. In the same way, we can understand the inter- 
face states at the junction between two TIs. We observe 
that the gap of the bulk band structure of the Hamil- 
tonian (1) closes at k ^ if the parameter A changes 
sign, and it closes at k = if the parameter B changes 
sign, which could be an indication of a topological phase 
transition. The origin of the interface states within the 
model can then also be understood from this closing of 
the bulk gap when A = at k ^ for case 1, and when 
B = at k = for case 2. When both parameters A and 
B change sign as in case 3, we end up with two TIs with 
equal helicity direction. Therefore, at first sight it might 
be surprising that interface states appear. However, from 
the previous discussion it is clear that the bulk gap closes 
at two different k values when both A and B change sign, 
which explains why the resulting spectrum in case 3 (Fig. 
4) is a combination of the spectra of case 1 (Fig. 2(a)) and 
case 2 (Fig. 3). This is illustrated with a schematic phase 
diagram in Fig. 5, where we show the four different sign 
combinations of A and B, and the (k, k z ) points at which 
the gap closes when the helicity changes orientation. 

Because of the invariance of the Hamiltonian (1) under 
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FIG. 6: (color online). Dispersion of the interface states for 
A (2) = -A w and B (2) = 3.0 eVA, where we have included 
the tachyonlike solutions; this corresponds to Fig. 3(a) of Ref. 
21. Eq. (4) is plotted for Til (red) and TI2 (green) together 
with the numerical results (dots). From Eq. (5) we find fc^ 1 ' m 
0.0853 and fc t (2) w 0.0381. 

any rotation about the z-axis, we find gapless interface 
states along any fc-direction. If we take into account 
terms up to order fc 3 , the full rotation symmetry is re- 
laxed to threefold C3 rotation symmetry of the Bi2Sc3 
crystal 20 , and we expect a gap to open in the interface 
spectrum along k x ^ or k y ^ directions, giving a total 
of six Dirac cones over the entire Brillouin zone. How- 
ever, the solution method becomes unpractical in this 
case, because the equation for A (A.l) contains odd pow- 
ers, and we need to find the roots of a depressed quartic 
equation. 

Now raises the question if materials can exist for which 
the parameters A and B changes sign with respect to 
their Bi2Sc3 values. In Ref. 20 it is argued that the sign 
of A is determined by the sign of the atomic SOC, while 
the sign of B is independent of the atomic SOC. How- 
ever, the atomic SOC parameter, given by the diago- 
nal matrix element of Hso = (2m 2 c 2 r) _1 (9f//9r)i • S, 
does not change sign for Bi 2 Se3-like TIs because the po- 
tential is always attractive for atoms 20 . Therefore, the 
helicity direction is always the same within the Bi2Sc3 
family of TIs. However, the model (1) is valid for a gen- 
eral TI, and therefore we believe that the interface be- 
tween TIs with opposite helicity is possible in principle, 
as was also mentioned in Ref. 19. Furthermore, there are 
other mechanisms besides SOC which can give rise to a 
band inversion. For example, strained bulk HgTc is a 3D 
TI 27 in which the band inversion is caused by other rel- 
ativistic corrections 28 . Also in typc-II InAs/GaSb het- 
crostructures, one can achieve a band inversion from 
confinement 29 . 



E. Tachyonlike interface states 

We found no signature of tachyons in this system as 
claimed in Ref. 21. The reason is the following. If the 
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A's used to construct the general solution (2) become 
four times degenerate, the solution is no longer valid be- 
cause it contains the same contribution twice. In this 
case, the matrix of the homogeneous system, defined by 
the boundary conditions (3), is always singular and one 
seemingly finds solutions. The extra degeneracy occurs 
at (E,k) points where the equation in A 2 has a repeated 
root (Eq. (A.l)). This happens when the discriminant 
vanishes which allowed us to obtain an expression for the 
'dispersion', 



E t {k) 



with 



r Cl 



Mn 
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± 



1 ^(C 2 -M 2 )f(k), 



2M( 



(4) 



This topological phase transition has been characterized 
by a topological invariant produced by the mirror sym- 
metry of the system in Rcf. 19. 

Also, we found that the tachyonlike dispersion pre- 
sented in Ref. 21 does not correspond to physical solu- 
tions. Rather, it is a consequence of implementing the 
wrong general solution at points on the (E, k) grid where 
the secular equation (A.l) has four repeated roots. This 
degeneracy occurs independent of boundary conditions 
and the same complication arises for a vacuum interface 
where one can analytically show that these solutions are 
unphysical. From the comparison of Eq. 4 with the nu- 
merical results in Fig. 6, and the fact that the tachyonic 
solutions vanish if we use the correct general solution at 
the degeneracy, we conclude that this system does not 
show any signature of tachyons. 



f{k) = B 4 + AB 2 M a M 1 + 4M 1 (B 2 M 2 - A 2 M x )k 2 
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k t = -r- \ 



B 2 + 4Af Mi 
Mi(A 2 Ah - B 2 M 2 



(5) 



In Fig. 6 we show the dispersion for the case which was 
also considered in Rcf. 21, = and B^ = 3.0 

cVA. We see that Eq. (4) exactly fits the tachyonlike 'dis- 
persions' obtained with a numerical calculation that does 
not take the degeneracy into account. In this case we get 
two tachyonlike 'dispersions' because the degeneracy oc- 
curs for A' 1 - 1 and A*- 2 -* for this parameter set. We know 
that the degeneracy occurs on a real (E, k) grid if k t is 
real. The same problem occurs if we consider a vacuum 
interface, because the degeneracy of A is independent of 
the boundary conditions. In this case, one can analyti- 
cally show that the tachyonic solutions are not physical. 

The correct general solution at the degeneracy is given 
by Eq. (A. 6) in the Appendix. Now we only find solu- 
tions at the points of intersection between the unphysical 
tachyonlike 'dispersion' and physical solutions. 



IV. CONCLUSION 



Appendix: General solution 

After the substitution k z — > —id z , the Schrodinger 
equation (H — E) <f) = is given by a coupled system of 
homogeneous second order differential equations that is 
solved with the ansatz 4>(z) ~ e Xz ^(E, k), where ^(E, k) 
is an eigenvector of the Hamiltonian. This yields a sys- 
tem of algebraic equations which have a nonzero solution 
if \H(k, —iX) — E\ = which we can write as the square 
of a biquadratic equation in A, given by 



L>iD 2 A 4 + (B 2 +D 1 {E- L 2 ) + D 2 (E - L x )) A 2 



(E — L±)(E - L 2 ) - A 2 k 2 =0, 



(A.l) 



with 



Li, 2 = C ± M + (C 2 ± M 2 ) k 2 , 
D h2 = C 1 ±M 1 . 

This gives four doubly degenerate X a (E,k) in general, 
and the eigenspace of X a is spanned by two eigenvectors 
^ s ,a{E, k) which can be chosen as 22 



In conclusion, we have studied interface states at the 
junction between two TIs with equally and oppositely 
oriented helicity direction. The origin of the gapless in- 
terface states can be understood from the closing of the 
bulk gap in the transition from one orientation of helic- 
ity to the opposite orientation at different wavevectors. 
Even if we relax the full rotation symmetry of the model, 
there should still be gapless states at k x = or k y = if 
we wish to avoid the scattering paradox of surface states 
at the junction between two TIs with opposite helicity. 
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Li+AA 2 
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Ak_ 
L 2 + D 2 X 2 D 

lBX a 



(A.2) 

The general solution is then given by Eq. (2). However, 
if the A are further degenerate, i.e. when Eq. (A.l) has a 
repeated root, there are only two distinct A that are both 
four times degenerate. In the last case, there are only two 
linear independent solutions for each A, given by e Az ^ / i 
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and e Xz ^>2- In order to find two more linear independent 
solutions, we first consider a general system of n homo- 
geneous second order ordinary differential equations: 

0" = Q0' + P<j), (A.3) 

where P and Q are constant n x n matrices. If we plug 
4>(z) = e Xz rj in Eq. (A.3), we find 

(A + XB- X 2 l) r] = 0, (A.4) 

which has a nonzero solution 77 ^ if and only if 
det (A + XB — A 2 /) = 0. In general, the solution of this 
equation is given by the roots of a polynomial of order 
2n in A. Now consider the case where there is a repeated 
root. If there are two linearly independent r) correspond- 
ing to the repeated root, there is no problem. When this 



is not the case, we need to look for another linear in- 
dependent solution. If we try 4>(z) = ze Xz r) + e Xz p, and 
equate terms of equal power in z, we find Eq. (A.4) again 
and additionally we have 

(A + XB - A 2 /) p = (2X1 - B) rj. (A.5) 

If we take the partial derivative of Eq. (A.4) with respect 
to A, and compare the result to Eq. (A.5), we find p = 
d\T]. If we use this result, we obtain the general solution 
when the A are four times degenerate: 

4>(z) = J2 eX " Z [ C :P*:f> + C ^+ 2 ( z + *:f>] ■ 

8,0=1,2 

(A.6) 
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